Transient resonances in the inspirals of point particles into black holes 
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We show that transient resonances occur in the two body problem in general relativity, for spinning 
black holes in close proximity to one another when one black hole is much more massive than the 
other. These resonances occur when the ratio of polar and radial orbital frequencies, which is slowly 
evolving under the influence of gravitational radiation reaction, passes through a low order rational 
number. At such points, the adiabatic approximation to the orbital evolution breaks down, and 
there is a brief but order unity correction to the inspiral rate. The resonances cause a perturbation 
to orbital phase of order a few tens of cycles for mass ratios ~ 10 -6 , make orbits more sensitive 
to changes in initial data (though not quite chaotic), and are genuine non-perturbative effects 
that are not seen at any order in a standard post-Newtonian expansion. Our results apply to an 
important potential source of gravitational waves, the gradual inspiral of white dwarfs, neutron 
stars, or black holes into much more massive black holes. Resonances effects will increase the 
computational challenge of accurately modeling these sources. 



Introduction: The dynamics of a two-body system emit- 
ting gravitational radiation is an important problem in 
general relativity. Binary systems of compact bodies 
undergo a radiation-reaction-driven inspiral until they 
merge. There are three different regimes in the parame- 
ter space of these systems: (i) The weak field, Newtonian 
regime r» 1, where r is the orbital separation in units 
where G = c = M = 1 and M is the total mass. This 
regime can be accurately modeled using post-Newtonian 
theory, which consists of an expansion in the small pa- 
rameter 1/r pp. (ii) The relativistic, equal mass regime 
r ~ 1, e ~ 1 (where e = n/M is the mass ratio with 
/i being the reduced mass), which must be treated using 
numerical relativity. Numerical relativists have recently 
succeeded for the first time in simulating binary black 
hole mergers [3j . (iii) The relativistic, extreme mass ratio 
regime r ~ 1, e <C 1, which is characterized by long, grad- 
ual inspirals on a timescale ~ £~ 1 , and for which compu- 
tational methods are currently under development. 

In this Letter, we show that in the relativistic, extreme 
mass ratio regime, there are qualitatively new aspects to 
the two-body problem in general relativity, namely the 
effects of transient resonances. While resonances are a 
common phenomenon in celestial mechanics when three 
or more objects are involved [3j, they can occur with just 
two objects in general relativity, due to its nonlinearity. 
They are a non-perturbative effect for highly relativistic 
sources and are not seen at any order in standard post- 
Newtonian expansions of inspiral solutions. Their exis- 
tence is closely related to the onset of chaotic dynamics, 
which has previously been shown to occur in general rel- 
ativity in other, cosmological contexts pt]. 

The resonances have direct observational relevance: 
Compact objects (1 < fx/M & < 10, where M Q is the 
Solar mass) inspiraling into much larger black holes are 
expected to be a key source for gravitational wave detec- 



tors. Advanced LIGO will potentially observe 3 — 30 such 
events per year, with 50 < M/Mq < 1000 [5], and future 
space-based detectors are expected to detect such inspi- 
rals with 10 4 < M/Mq < 10 7 out to cosmological dis- 
tances at a rate of ~ 50 per year [5] . The observed grav- 
itational wave signal will be rich in information. For ex- 
ample, one will be able to extract a map of the spacetime 
geometry of the central object and test if it is consistent 
with general relativity's predictions for a black hole [5J[7]- 
Such mapping will require accurate theoretical models of 
the gravitational waveforms, which remain phase coher- 
ent with the true waveforms to an accuracy of ~ 1 cycle 
over the large number ~ e -1 ~ Mj fj, ~ 10 2 — 10 6 of cycles 
of inspiral. Over the past decade there has been a sig- 
nificant research effort aimed at providing such accurate 
models [8 . Resonances will complicate this enterprise, 
as we discuss below. 

Method of Analysis: Over timescales short compared to 
the dephasing time ~ e^ 1 / 2 , inspirals can be accurately 
modeled using black hole perturbation theory, with e as 
the expansion parameter. The leading order motion is 
geodesic motion on the background Kerr metric. At the 
next order the motion is corrected by the particle's self- 
force or radiation reaction force, for which a formal ex- 
pression is known [5], and which has been computed ex- 
plicitly in special cases; see, e.g., the review |S]. Over 
the longer inspiral timescale ~ e^ 1 , it is necessary to 
augment these methods with two-timescale expansions 
which are currently under development [101 lllj . In this 
framework the leading order motion is an adiabatic in- 
spiral, and there are various post-adiabatic corrections. 

Geodesic motion in the Kerr spacetime is an integrable 
dynamical system, and it is useful to use the the corre- 
sponding generalized action-angle variables to parame- 
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terize the inspiral. The resulting equations are [T7] : 
^ = o Ja (3)+egi 1 H qe ,q r ,J) + 0(e 2 ), (la) 
^ = eGl 1 \ qe , qr ,J)+s 2 G^( qe ,qr,J) + 0(e 3 ). (lb) 

Here r is Mino time [13] and J„ are the con- 
served integrals of geodesic motion given by J„ = 
(E I fx, L z /fx, Q/fJ, . ), where -E is the energy, L z is the angu- 
lar momentum, and Q the Carter constant. The variables 
la = (<Zti <7r, Qe, q <j)) are generalized angle variables conju- 
gate to Mino time 10J. The right hand sides at O(e ) 
describe geodesic motion, with fundamental frequencies 
oj r , tug and uj^. The forcing functions g„ , G„ and G^ 
are due to the first order and second order self-forces, and 

(1 2) 

are 27r-periodic in qg and q r . The piece of G v : that is 
even under qg — > 2ir — qg, q r 2ir — q r , and the piece 
of g a X ^ that is odd, are the dissipative self-force, and the 
remaining piece is the conservative self-force [10] . 

In the limit e — > 0, the leading order solutions to Eqs. 
([I]) are given by the following adiabatic prescription |10j : 
Drop the forcing terms and G^ , and replace G^ by 
its average over the 2-torus parameterized by q g, qr . It is 
now known how to evaluate this averaged force explicitly 
[12] [13] , although generic adiabatic inspirals have not yet 
been computed numerically. 

Consider now post-adiabatic effects. The dynamical 
system ([I]) consists of a perturbed, integrable Hamilto- 
nian system. Resonances in this general type of system 
have been studied in detail and are well understood [H] , 
and we can apply the general theory to the present con- 
text. The existence of resonances in this context has pre- 
viously been suggested by Refs. [TSJIIII]- We will present 
three different treatments of the resonances: (i) An intu- 
itive, order of magnitude discussion, which is sufficient to 
deduce their key properties; (ii) A numerical treatment; 
and (iii) A sketch of a formal analytic derivation. A more 
detailed treatment will be presented in Ref. [17] , 

Order of Magnitude Estimates: Suppose that we have an 
adiabatic solution, which will be of the form q a (r, e) — 
ip a (sT)/e, J„(t,s) = J u {£t). Consider now the post- 
adiabatic correction terms in Eqs. ([I]), near some ar- 
bitrarily chosen point r = 0. We expand q g as q g = 
qgo + LjgoT + dignT 2 + 0(t 3 ), where subscripts denote 
evaluations at r = 0, and we expand q r similarly. We also 
expand G^ as a double Fourier series: Gv{qe, qr, J) = 

Efe.n G i 1) kJ J ) et{kqe+nqr ^ where the 00 term is the ana- 
batic approximation, and the remaining terms drive post- 
adiabatic effects. Inserting the expansions of qg and q r , 
we find for the phase of the (fc, n) Fourier component 



(constant) + (kojgo + noj r o)T + (kihgo + n6j r n)T 2 



(2) 



Normally, the second term is nonzero and thus the force 
oscillates on a timescale ~ 1, much shorter than the in- 
spiral timescale ~ 1/e, and so the force averages to zero. 



However, when the resonance condition kcugo + nuj r o = 
is satisfied, the (k, n) force is slowly varying and can- 
not be neglected, and so gives an order-unity correc- 
tion to the right hand side of Eq. (lb). The duration 
of the resonance is given by the third term in |2| to be 
r rcs ~ \j\fvCj ~ 1/^/ve, where v = |fc| + |n| is the order of 
the resonance; after times longer than this the quadratic 
term causes the force to oscillate and again average to 
zero. The net change in the action variables J v is there- 



this change causes a phase error A(f> that accumulates 
over an inspiral, of order the total inspiral phase ~ 1/e 
times A J/ J ~ which gives A</> ~ 1/^/ve. 

This discussion allows us to deduce several key prop- 
erties of the resonances. First, corrections to the gravita- 
tional wave signal's phase due to resonance effects scale 
as the square root of the inverse of mass of the small body. 
These corrections thus become large in the extreme-mass- 
ratio limit, dominating over all other post-adiabatic ef- 
fects, which scale as e° ~ 1. 

Second, they occur when ui r /ojg is a low order rational 
number. There is a simple geometric picture correspond- 
ing to this condition [TBI HH] : the geodesic orbits do not 
ergodically fill out the [qg,q r ) torus in space as generic 
geodesies orbits do, but instead form a 1 dimensional 
curve on the torus. This implies that the time-averaged 
forces for these orbits are not given by an average over 
the torus, unlike the case for generic orbits. 

Third, they occur only for non-circular, non-equatorial 
orbits about spinning black holes. For other cases, the 
forcing terms Gu depend only on q g, or only on g r , but 
not both together, and thus the Fourier coefficient G„ kn 
will vanish for any resonance. 

Fourth, they are driven only by the spin-dependent 
part of the self-force, for the same reason: spherical sym- 
metry forbids a dependence on q g in the zero-spin limit. 

Fifth, they appear to be driven only by the dissipa- 
tive part of the self-force, and not by the conservative 
part, again because the forcing terms do not depend on 
both q g and g r . We have verified that this is the case up 
to the post-Newtonian order that spin-dependent terms 
have been computed [19j . and we conjecture that it is true 
to all orders. The reason that this occurs is that the con- 
servative sector of post-Newtonian theory admits three 
independent conserved angular momentum components; 
the ambiguities in the definition of angular momentum 
are associated with radiation, in the dissipative sector. 
As a consequence, the perturbed conservative motion is 
integrable to leading order in s, and an integrable per- 
turbation to a Hamiltonian cannot drive resonances. 

Sixth, although the resonance is directly driven only 
by dissipative, spin dependent self-force, computing res- 
onance effects requires the conservative piece of the first 
order self-force and the averaged, dissipative piece of the 
second order self-force. Those pieces will cause O(l) cor- 
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rections to the phases over a complete inspiral [TU], and 
the kicks A produced during the resonance depend on 
the O(l) phases at the start of the resonance. 

Seventh, resonances give rise to increased sensitivity 
to initial conditions, analogous to chaos but not as ex- 
treme as chaos, because at a resonance information flows 
from a higher to a lower order in the perturbation ex- 
pansion. For example, we have argued that changes to 
the phases at 0(1) prior to the resonance will affect the 
post-resonance phasing at 0(l/y/e). Similarly changes to 
the phases at 0(^/e) before resonance will produce 0(1) 
changes afterwards. With several successive resonances, 
a sensitive dependence on initial conditions could arise. 

Numerical Integrations: The scaling relation A</> oc 1/ yfe 
suggests the possibility of phase errors large compared to 
unity that impede the detection of the gravitational wave 
signal. To investigate this possibility, we numerically in- 
tegrated the exact Kerr geodesic equations supplemented 
with approximate post-Newtonian forcing terms. While 
several such approximate inspirals have been computed 
previously [20] . none have encountered resonances, be- 
cause resonances require non-circular, non-equatorial or- 
bits about a spinning black hole with non orbit-averaged 
forces, which have not been simulated before. 

For the numerical integrations we use instead of q a the 
variables q a = (q tl q r , qe,q<j,) = (t, ip, \, 4>) where if> and x 
are the angular variables for r and 9 motion defined in 
Ref. [33]. The equations of motion in these variables 
are 

t, T = u t (q e ,q r ,J), <p T = cd^>(q$, q r , J), (3a) 
qe,r = tie(q~0, J) + shg(qg,q r , J) + 0(e 2 ), (3b) 
<Zr,r = w r (g r ,J)+e^ 1) (q e ,g r ,J) + 0(e 2 ), (3c) 
J VtT = eHi 1 \q e ,q r ,3)+0(e 2 ). (3d) 

Here r is Mino time [13j . the frequencies Co are given in 
[15] . and ha and are given in terms of the compo- 
nents of the 4-acceleration in |21j . 

We parameterize the three independent components of 
the acceleration in the following way: a a = a r e < ?+a e e? + 

a.±e a /3lS u l3 eJe 5 -+(a r Uf+a e Ug)u a , where u is the 4- velocity 
and Cf and eg are unit vectors in the directions of d r and 
dg. We compute the dissipative pieces of a r , a e and a± 
from the results of 22J, as functions of f = r + a 2 /(4r), 
E n = E — 1, and K = Q + a 2 : L 2 + a 2 E n , and then expand 
to 0(a 2 ) and to the leading post-Newtonian order at each 
order in a [17] ■ We also add the conservative component, 
expressed similarly and computed to 0(a) and to the 
leading post-Newtonian order, taken from Refs. [55J; see 
Ref. [T7] for details. 

We numerically integrate Eqs. ([3]) twice, once using 
the adiabatic prescription, and once exactly, and then 
subtract at fixed t to obtain the post-adiabatic effects. 
The adiabatic prescription involves numerically integrat- 
ing the right hand sides over the torus parameterized 
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FIG. 1: [Top] The adiabatic inspiral computed from our ap- 
proximate post-Newtonian self-force, for a mass ratio e = 
H/M = 3 x 10~ 6 , with black hole spin parameter a = 0.95, 
with initial conditions semilatus rectum p = 9.0M, eccentric- 
ity e = 0.7, and orbital inclination 9- mc ~ 1.20. The bottom 
curve is e, the middle curve is #i nc , and the top curve is ratio of 
frequencies wg/ccv, shown as functions of p. [Middle] The fluc- 
tuating, dissipative part of the first order self-force causes a 
strong resonance when u)e/uj r — 3/2 at p — 8.495. Shown are 
the corrections to the energy E, angular momentum L z and 
Carter constant Q, as functions of p, scaled to the ir values at 
resonance, and divided by the square root yj Jl/M of the mass 
ratio. The sudden jumps at the resonance are apparent, with 
the largest occurring for the Carter constant. [Bottom] The 
lower curve is the correction to the number of cycles <j>/(2ir) 
of azimuthal phase of the inspiral caused by the fluctuating, 
dissipative part of the first order self-force. The sharp down- 
ward kick due to the resonance at p = 8.495 can be clearly 
seen. The resonant corrections to the number of cycles of r 
and 9 motion are similar. These phase shifts scale as \/M//x. 
The upper curve is the post-adiabatic phase correction due 
to the conservative piece of the first order self-force, which is 
considerably smaller and is independent of the mass ratio. 



by qg,q r at each time step, where q r = F r (q r ) / ' F r (2ir) , 
F r (q r ) — J 9r dq r /[ui r (q r , J)], with a similar formula for 
qg . This is numerically time consuming, but the adiabatic 
integration can take timesteps on the inspiral timescale 
~ 1/s rather than the dynamical timescale ~ 1. 

Typical results are shown in Fig. [l] which shows the 
adiabatic inspiral for a mass ratio of e — 3x 10~ 6 with a = 
0.95, in terms of the relativistic eccentricity e, semilatus 
rectum p and orbital inclination f? inc , which are functions 
of E, L z and Q [53]. This example has a strong resonance 
at LOg/uj r — 3/2, that generates jumps in the conserved 
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quantities of order a few percent times yfe, and causes 
phase errors over the inspiral of order 20 cycles. Phase 
errors of this magnitude will be a significant impediment 
to signal detection with matched filtering. We find that 
the resonance effects are dominated by the 0(a 2 ) terms, 
and the effect of the O(a) terms are small. 

Analytic Derivation: In terms of the slow time variable 
f = st, the solutions of the dynamical system ([I]) away 
from resonances can be expressed as an asymptotic ex- 
pansion in e at fixed f [TUl [H] : 



C L 

J v {r,e) = J^(f) + V^jy /2) (f) + 0(e). 



. (4a) 
(4b) 



The leading order terms give the adiabatic approximation 



described above, and satisfy [TU] ip a \ = lu c 

(G[z 1 ' ) )[J7' < ' '], where the angular brackets denote an aver- 
age over the (q r ,qe) torus. The subleading, post-1/2- 
adiabatic order terms satisfy J^ 1 / 2 ^ — {Gv = 

AJ^ 1/2) <5(t), ip£t 2) = u a j^jtl 1/2 \ where the (5-function 
source term arises at a resonance, taken to occur at r = 0. 

Near the resonance we use an ansatz for the solutions 
which is an asymptotic expansion in ^fe at fixed f = \/et, 
and then match these solutions onto pre-resonance and 
post-resonance solutions of the form Q [H] . To linear 
order in the force Fourier coefficients (an approximation 
which is valid here to within a few percent |17|). the 
jumps in the action variables for a resonance (fc, n) can 
be computed by substituting the adiabatic solutions into 
the right hand side of Eqs. ([I]) and solving for the per- 
turbation to the action variables. The result is 



a #n = £ . 



2tt 



as 



exp 



sgn(as) — +«sxrc 



G 



(i) 



where x rcs = kqg + nq r , a = kuig.f + nw r .f and all quan- 
tities are evaluated at the resonance r = using the 
adiabatic solution. In Ref. [T7] we give the exact ex- 
pression for this quantity that does not linearize in the 
force Fourier coefficients. We note that evaluating the 
phase Xres requires knowledge of the second subleading, 
0(1) phase in Eq. (4a), which in turn requires knowl- 



edge of the force components g„ \ G^ and (G'if 1 ) in Eq. 
([T]) [TU]. In addition, to obtain the phase to 0(1) accu- 
racy after the resonance, it is necessary to also compute 
the subleading, O(e) jumps in J„ and 0(1) jumps in q ai 
which are given in [17] , 

Discussion: The dynamics of binary systems in general 
relativity is richer than had been appreciated. Transient 
resonances occurring during the inspiral invalidate the 
adiabatic approximation and give rise to corrections to 
the orbital phase that can be large compared to unity. It 
will be necessary to incorporate resonances into theoreti- 
cal models of the gravitational waveforms for inspirals of 



■y(2)\ 



compact objects into massive black holes, an important 
gravitational wave source. This will require knowledge 
of the second order gravitational self force and will be 
challenging. 
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